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ABSTRACT : 

This report deals with the aerodynamic design of an axial compressor 
stage with symmetrical bladings for a research program to investigate 
tip clearance effects in the three-stage compressor of the Turbo- 
Propulsion Laboratory, NPS. It establishes the blading data and the 
stage performance with an iterative three-dimensional approach, and 
gives design criteria for the drive and the flow measuring device of 
the test unit. The calculated distributions of the flow properties in 
the stage will be used for future comparisons with test data. 
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1 . INTRODUCTION 

The three-stage axial flow compressor test rig of the Turbo- 
Propulsion Laboratory, Department of Aeronautics, NPS, will be used to 
carry out a Navy sponsored research program on three-dimensional flow 
phenomena in compressor bladings. This compressor with a tip diameter 
of 36 in. and a hub/tip ratio of 0.6 is presently equipped with lightly 
loaded bladings which are of the free-vortex type. A 50 HP two-speed 
electric motor is now directly coupled to the compressor. Further 
details of the test rig are given in Ref. 1 (See Bibliography at end of 
text) . 

During the first phase of the research program, so-called symmetrical 
bladings with high aerodynamic loading will be tested. Further, in 
order to increase the pressure rise in the compressor, primarily to 
improve the measuring accuracy, the present motor will be replaced by a 
150 HP motor that was obtained from surplus. This motor has the follow- 
ing characteristics: 

Make: General Dynamics/Electro Dynamics, Bayonne, N. J. 

Frame: 505 Y - Type TN 

Serial No.: 70105116 A6 

Cont. Load: 150 HP, 1180 rpm, 4A0V/3 phase, 190 Amps 

It is planned also to make tests with one, two, or three stages, 
preferable in such a manner that the full power of the motor can be 
utilized in all three cases, requiring operations at different speeds. 

It is intended to make this possible by arranging a V-belt drive with 
different replaceable pulleys. 

The objective of this study is to establish the aerodynamic design 
of the new blading, and to determine the maximum speeds at which these 
stages can operate at a driving power of 150 HP. These data will then 
be used in a future report to design the V-belt drive and the nozzle for 
the measuring of the compressor flow rate, since Ref. 1 shows that the 
presently used measuring method gives doubtful results. 

The blade loadings will be such that the NASA diffusion factor does 
not exceed a value of 0.5 for the rotor, or 0.6 for the stator. All 
three stages will be identical. The flow conditions ahead of and after 
the blade rows will be determined with the methods of Ref. 2. 
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2. FIRST APPROXIMATION OF THREE-DIMENSIONAL FLOW PATTERN. 

For stages with equal energy input from hub to tip without imposed 

radial energy gradients, the change of the axial velocity component in 

radial direction can be obtained from Eq. 16(50) of Ref. 2. For a first 

approximation it is assumed that the axial components V ^ and V j ahead 

of and after the rotor blades are equal at the mean radius R , and that 

m 

the radial shift of the stream surfaces can be ignored. The resulting 

velocity diagram of the blading at the arithmetic mean radius R is 

m 

shown in Fig. 1. Figure 2 represents the conditions for a cylindrical 
stream surface at an arbitrary radius R between hub and tip. The actual 
velocity diagram for a stream surface that has the radii R^ and R^ ahead 
of and after the rotor, respectively , is given by Fig. 3. The symbols 
used in these figures, and in the following derivations, are defined in 
Ref. 2. 

From Eq. 16 (42) of Ref. 2 and by Fig. 1 



K-RAV-R AV « R V (tang, - tang, ) 
u m urn m m lm 2m 



Also 



wR - W , + V , « W , +W 0 = V (tang, + tang. ) 

m ulm ulm ulm u2m m lm 2m 



and 



u R 



m 



m tang, + tang„ 
lm zm 



thus 

K - w R 2 X 
m 

where 



X - tan6 lm - tane 2m 
tanS lm + CanS 2a, 
From Eq. 16 (52) of Ref. 2 



V ul . _ fm X 1 R 

b)R R 2 2 R 

m m 
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m 111 
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From Eq. 16 (53) of Ref. 2 
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m 2 2 
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Introduced into Eq. 16 (50) of Ref. 2 for 
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From Eqs. 8, 7 and 2 
2 
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where 



(tan6 lm + tan6 2m >' 



( 11 ) 



From Eqs. 9, 7, 2, and 10 
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Thus, the axial velocity components V ^ and V £ at different radius 

ratios R/R can be determined with Eqs. 11 and 12 if the angles £5, 
ra lm 



and are specified at the mean radius R 

' 1 



m 



From Fig. 2 
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ul 



al 
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With Eq. 5 
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Similarly, with Eq. 6 
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Also from Fig. 2 and Eq. 5 
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and with Eq. 6 



tan a. 



V 



u2 
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R , R 
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With these flow angles the relative and absolute velocities W^, V^, 

and V , V are known as multiples of coR . 
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3. PERMISSIBLE FLOW DEFLECTIONS 



The NASA diffusion factor D r of the rotor flow is defined by 



[Ref. 6] 
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W AW 

d r " 1 ' ^ + ‘ 15 ?T 



( 17 ) 



where o R is the rotor solidity, 
stator flow is 

V- AV 

° s " 1 ‘ \ + T ^r r i 



The diffusion factor D of the 

s 



(18) 



Assuming that the stator solidity o equals o_, and that the chord 

S K 

of the stator and rotor blades does not change along the radius, 
there is 



m 



o_ * o * o _ 
R s m R 



(19) 



where o is the solidity at R . 
m J m 

From Eqs. 1 and 3 



AV 



u 



AW - K 
U R 



R R m X 

“ R mlT 



( 20 ) 



Then from Eq. 17 with 19 



V -/uiR cosB l X cos Ei 
n = 1 a2 m ■ 1 

R “ V_,/o)R_ co*6~ 2 a (V . /to R ) 



al m 



m al m 



( 21 ) 



Also by Eqs. 18, 19, and 20 



V ,/wR cosa. X cosa- 

n ml— fl l m 4. z 

s V _/wR cosa. 2 o (V 0 /uR ) 

a2 m 1 m a2 m 



( 22 ) 



For the chosen conditions the diffusion factors D n and D at the mean 

R s 

radius R m are equal, say, equal to D , where by Eqs. 21, 2, and 4 



cos8 

D « 1 

m cosB 



lm + i tan6 l„ - tan6 2» ) coa6 lm 



2m 



2 o 



(23) 



m 
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Equation 23 can be rearranged to calculate directly the flow angle 
& 2 m f° r specified values of B^ ffi , D m > and o^. Equation 23 rewritten 
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tanB. tanB 0 

lm zm 
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lm cos g 
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m m 



lm 



then 
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2m 
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m (-) V 
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(24) 



(25) 



It is possible also to determine the deflection angle 



Ag 

from Eq. 
2o 

m 



* S lm " B 2m 

23. This relation rearranged 

cosBo (1 - D ) = 2o cosB. - (tanB. - tanB,, )cosB, cosB 

zm m m lm lm zm lm 
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m lm lm 2m lm 2m 



(26) 



2m 
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m lm lm 2m / 
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m lm 
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Since from Eq. 26 



6 2m - 6 i m - Ae 



there is 
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m in im lm m ±m 
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there is 
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m m lm 
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mm lm J m 
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(28) 



For chosen values B, » a , and D the flow angles B~ can therefore 

lm m m Zm 

be determined. The angles B^ m and B2 m then establish the quantities 
X and Y of Eqs. 4 and 11, respectively. Introduced into Eqs. 10 and 
12 these values can be used to obtain the changes of the axial com- 
ponents and in radial direction. With the angles B^ and 
from Eqs. 13 and 14 it is then possible to calculate the diffusion 
factors D and D at particular radii with Eqs. 21 and 22. 

R S 

4. APPROXIMATE OPERATING CONDITIONS OF COMPRESSOR. 

The following relations will be used to establish the main 
parameters of the blading, which later on will be investigated in more 
detail. 

If it is assumed that V is equal to the average through-flow 

m 

velocity in the compressor annulus, the volume flow rate Q is 

Q = A V 
cm 

where A is the cross-sectional area of the annulus. With Eq. 2 
c 



Q * A wR 
cm 



tan6 lm + tan6 2m 



[ft /sec] 



(29) 
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The driving moment M for one stage is 



M “ mR._ aV » mK 

fn **> 



[ft-ib] 



Because of the small pressure rise in the compressor the mass flow 

O 

rate m can be taken as 
m - Q p 

where p is the air density at ambient conditions. For z stages 
(z ■ 1, 2, or 3) the necessary horsepower HP is 



HP 



u M_ 

550 



3 3 

Po co R A (tanB, -tanB 0 ) 
me lm 2 m 

550 



At 14.7 psia and 60 F 

p - 2.371 (10* 3 ) [slug/ft 3 ] 

O 

With the rotational speed N in rpm 

it N 
“ ' 3CT 

For a tip diameter of 3 feet and a hub/ tip ratio of 0.6 

A - 7 - 3 2 (l-0.6 2 ) - 4.524 ft 2 
c 4 



R - 1.5 — + -°— I = 1.2 ft 
m 2 

Then, from Eq. 29 



Q “ 



(0.568) N 



tanB. + tanB 0 
lm 2 m 



[ft /sec] 



From Eq. 30 



HP 



3 

<38 - 7)z (io5o) <“'l. ■ “V 



Assuming that 135 HP are absorbed at the design point, to have 
available 150 HP at overload conditions, the design operating speed 
is obtained from 



N 



1516.6 



z (tanB lm - tanB 2m ) 



1 * 



[rpm] 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 
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If AP fc is the rise of the total pressure in the stage, and 
the total-to-total stage efficiency, there is 



Q AP 

=• Mco - HP (550) 

T l 



With Eqs. 30, 31 and 3 



( 36 ) 



AP t = Mu> * p Ku « p oj 2 R 2 X (37) 

T~ • m 

ou 

For^driving power of 135 HP, the pressure rise AP^ is obtained 
directly from the second equality of Eq . 36. With Eq . 33, 



or 



AP tang. + tang 

-r ‘ 130 - 72 mm !lb/£t 1 



AP. 



25.14 



tang, + tang„ 
im zm 



r| (N/1000) 

5. DETERMINATION OF BLADING PARAMETERS. 



[in.H 2 0] 



(38) 



Most desirable for the present tests are bladings that produce 
the highest pressure rise at the lowest speed of rotation. Such 
bladings will also have low flow rates. Figures 4 and 5 show that the 
angles g^ should be as high as possible for these conditions. 

Figure 6 indicates that the speed of rotation is not changing greatly 
with g^ and that its magnitude is primarily a function of the dif- 
fusion factor D^, which should be as high as possible. 

Figure 7 shows the flow deflections Ag * g^ - for different 

angles g.. and diffusion factors D at a solidity a “1, which 
lm tn m 

establish the quantities X and Y of Eqs. 4 and 11, respectively. These 
values introduced into Eq . 12 show that the axial velocity 
becomes imaginary for g^ greater than a particular angle which is 

depending on D . These limits of g- * for V 9 * 0 are: 
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D = 0.30 0.35 0.40 0.45 

m 

B lm * - 45.7° 47.2° 47.8° 48.4° 

Hence the flow angles 0^ cannot exceed about 45°. 

Table I gives the flow conditions along the radius for 0- * 45° 

lm 

and « 0.4 at the mean radius R , as calculated by the previously 
established relations* It can be noted that the diffusion factor D 

R 

at the rotor tip is in excess of 0.56, hence unacceptably high. For 

0^ m * 45° and * 0.35 the resulting data are listed in Table II. 

Although D at the rotor tip is about 0.5, the large deceleration of 

the axial velocity components at the outer radius from » 0.4445 

to ■ 0. 2629 f together with the negative flow deflection A8 at 

this station,make this blading undesirable for the present purposes. 

The angle 8^ m is therefore taken as 40° and Table III lists the data 

obtained for D *0.4. For this value the diffusion factor D_ at the 
m R 

rotor tip is larger than 0.5. Repeating the calculations for 8- = 40° 

lm 

and * 0.35 gave the flow conditions of Table IV which seem acceptable 

in all respects. These parameters are therefore chosen for a more 

detailed investigation of the flow conditions by taking account of the 

radial shift of the stream surfaces from rotor inlet to rotor discharge 

with the resulting velocity diagram of Fig. 3. 

6. METHOD FOR BETTER APPROXIMATION OF FLOW CONDITIONS. 

The quantities listed in Table IV will be called the first 

approximation of the flow conditions. Figure 8 shows the velocity 

ratios V ,/(a)R ) and V 0 /(u)R ) of this approximation as functions of the 
al m az m 

radius ratio R/R^. These values were obtained by assuming that and 
a t R m are equal, and identical with the quantity of Eqs. 8 and 9. 
Moreover a fluid particle entering the rotor at a radius R^ is supposed 
to leave the rotor at the same radius, or R2 ■ R^ ■ R. This condition 
led to Eq. 1 which formulates that all fluid particles receive the same 
energy increase AH in the rotor. However if R^ and Rj differ, this 
requirement must be expressed by 

AH * u)(R_ V _ - R, V ) * Constant (39) 

2 u2 1 ul 
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If it is assumed that the energy increase AH in the stage is 
that of the first approximation, there is from Eqs. 1 and 3 



AH - uK - go 2 R 2 X (40) 

m 

where, from Table IV, 

X « 0.36578 



Now, if the tangential components along the radius R^ ahead of the 
rotor are taken to be those of the first approximation, namely, by Eq. 5 




u)R 

m 



J1 X 1 ^1 

R. 2 2 R 

1 m 



(Al) 



the distribution of V ^ along the radius after the rotor is from 
Eqs. 39, AO and Al, 




(A2) 



Thus, to satisfy Eq. 39, the components V £ can be determined only if 

the corresponding radii and R^ of the stream surfaces are known. 

The results of the first approximation will be used to establish an 

approximate relationship between R^ and R^ by applying the equation of 

continuity. In particular, for a chosen value of R^ the corresponding 

radius R^ is found from the condition that the same mass flow rate 

has to occur between R 0 and the hub radius R 0 , , as between the radius 

2. 2n 

and the hub radius T ^e ac *°P tec * procedure is explained in the 

next paragraph. With R^ known as a function of R^, the distribution of 
the components V ^ in radial direction are obtained from Eq. A2. The 
axial components along R^ can then be determined by solving Eq. 16 

(Al) of Ref. 2, or 



3 (V 

Tr 




Z 



2 



2 V 



u2 



j Vul 

3R„ 



(43) 
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if the curvatures of the generatrices of the stream surfaces and radial 
entropy gradients are ignored. The data obtained in this manner from 
Eq. 43 represent the so-called second approximation. Applied to the 
equation of continuity they permit a check of the relationship be- 
tween and R 2 * If this interdependence does not correspond to that 
obtained from the first approximation, it will be necessary to carry 
out further iterations. 

Evidently, if curvature effects and entropy gradients are ignored, 
the distribution of V ^ along for V ^ of Eq. 41 is not affected by 



the radial shift 


from R- and 


R ? , hence the 


flow conditions ahead of 


the rotor are those given in 


Table 


IV. 




The relations for the ; 


second 


and , if 


necessary, further approxi 


mations, can be 


simplified by introducing the following dimensionless 


quantities : 


R„ 






R_ 




1 






2 


r l 


R 




r 2 * 


R 




m 






m 


* 


V al 




* 




V al 


(dR 




V a2 " 


uR 




m 






m 


* 






* 

IT . 


, !u2 


V ul 


uR 


* 


V u2 


wR 




m 






m 


* 

\r 


, JSi. 






jm2 


ml 


wR 

m 


i 


V m2 


wR 

m 



Then, Eq. 42 rewritten is 



u2 2r 



h ( x + r i 2 ) 



(45) 



and Eq. 43 becomes 



8 ( V a2 > 
3r„ 



S ♦ (?) 

r 2 V r 2 



x ) 

2 Tr~ 



Since X is a constant. 
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»( V a2* ) 



3r, 



X + r. 



fl 

3r, 



- f(r 2 ) (46) 



Hence, if is knovra as a function of r^, the function f(r 2 ) can be 
established. Denoting the as yet unknown value of V^ 2 at 1 1 by 
V . , the integration of Eq. 46 gives 



( V a 2 ‘) 2 ’ K 2 *) 2 -/ f(r 2> dr 2 < 47 

l 

A 

The quantity must be obtained by iterations with the help 

of the equation of continuity. Figure 4 shows that the pressure rise 
per stage does not exceed 16 inches of water. In a stage with 50 
percent reaction the pressure rise in the rotor is then about 8 inches 
of water, giving a pressure ratio of about 1.02 for an inlet pressure 
of 14.7 psia. The ratio $2^1 t ^ xe mass density a fter and ahead of 
the rotor will therefore not exceed a value of 1.014. 

The volume flow rate ahead of the rotor between the hub 
radius and an arbitrary radius is 

«1 - 2 "/\ V al dR l 
R lh 

With the quantities of Eq . 44 let 



Q i “ 2 ™ V 5 Q i* 

by introducing the dimensionless quantity 

r. 



(48) 



/ 



r i va i dr i 



(49) 



R lh /R m 



The total volume flow rate Q. ahead of the rotor Is 

lmax 

3 * 

Q, ■ 2iru> R Q, 
lmax m lmax 



(50) 
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where Q^ max is obtained from Eq. 49 for the upper limit r^ * R^/R^ 
of the integral, where R is the tip radius. 

Similarly, for the conditions after the rotor the volume flow 



rate Q 2 between R^ and R 2 is 
Q 2 = 2™R m 3 Q 2 * 



( 51 ) 



with Q, 



/ r 2 l 



42 * dr 2 



(52) 



R 2h /K m 



The total volume flow rate 0. at station 2 is 

^2raax 



^2max 



3 * 

2tto) R Q 0 

m 2max 



(53) 



The quantity Q 2 max is obtained from Eq . 52 for the upper limit 

r 0 * R 0 /R °f the integral. 

L It m 

The equality of the mass flow rates at stations (1) and (2) 
can be expressed by 



P, Q * p, Q 0 k ’ 
1 max 2 2max b 



where p^, p 2 are the average mass densities at (1) and (2), and kg’ 

is the so-called blockage factor that takes account of the increased 
regions of reduced velocity in the so-called wall boundary layers 
near the hub and the tip at station (2), compared to those at station 
(1). Then, with Eqs. 50 and 53 



Umax 



Umax 



Umax 



(p 2 /p l )k B ’ 



(54) 



Experience shows that the blockage factor k. ’ varies between 0.96 and 

” * 

0.97. Hence, since p p ¥ 1.014, the quantity Q 2 max be 

determined by Eq. 54 with^a factor k * (1 .014) (0.965) ■ 0.98 for the 
* ® * 
value of Q, obtained from Eq. 49. This value of 0- is used to 
Umax n 2max 
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X + r. 



^1 

3r, 



- f(r 2 ) 



(46) 



Hence, if r.. is known as a function of r_, the function f(r») can be 
1 1 * L 
established. Denoting the as yet unknown value of V , at r_ * 1 by 
* a2 2 

V . , the integration of Eq. 46 gives 



( V a 2 ‘) 2 • K 2 V-/ f < r 2> dr 2 < 47 

1 

£ 

The quantity V j must be obtained by iterations with the help 
of the equation of continuity. Figure 4 shows that the pressure rise 
per stage does not exceed 16 inches of water. In a stage with 50 
percent reaction the pressure rise in the rotor is then about 8 inches 
of water, giving a pressure ratio of about 1.02 for an inlet pressure 
of 14.7 psia. The ratio t ^ ie mass density after and ahead of 

the rotor will therefore not exceed a value of 1.014. 

The volume flow rate ahead of the rotor between the hub 
radius and an arbitrary radius is 

«1 ■ 2 ’) lR l V al dR l 
R lh 

With the quantities of Eq. 44 let 



3 * 

Q, ■ 2 ttw R Q, 

X 1 m 1 

by introducing the dimensionless quantity 



(48) 



/ 



r i Va i dr i 



(49) 



R lh /E „ 



The total volume flow rate Q, ahead of the rotor is 

lmax 



lmax 



3 * 

2ttu) R Q, 

m lmax 



(50) 



13 



where Q^ max is obtained from Eq. 49 for the upper limit r^ * R^ /R 

of the integral, where R^is the tip radius. 

Similarly for the conditions after the rotor the volume flow 

rate Q 0 between R 0 , and R 0 is 
L In l 

Q 2 = 27ru)R m 3 Q 2 * (51) 



with Q, 



/. r > l 



k2 dr 2 



(52) 



R 2b /E „ 



The total volume flow rate at station 2 is 

Umax 



2max 



2ttuj R 



m 



Umax 



(53) 



* 

The quantity Q 2 max obtained from Eq . 52 for the upper limit 

r„ * R_ /R of the integral. 

L it m 

The equality of the mass flow rates at stations (1) and (2) 
can be expressed by 



P i Q 
1 max 



p 2 



2max 



k 



b 



where p^, are the average mass densities at (1) and (2), and kg* 



is the so-called blockage factor that takes account of the increased 
regions of reduced velocity in the so-called wall boundary layers 
near the hub and the tip at station (2), compared to those at station 
(1) . Then, with Eqs. 50 and 53 



Umax 



Umax 



( p 2 /p 1 )I'b’ 



Umax 



(54) 



Experience shows that the blockage factor k. * varies between 0.96 and 

— _ B * 

0.97. Hence, since p^/ P =* 1.014, the quantity Q 2 max w iH be 

determined by Eq . 54 with^a factor k * (1.0i4) (0.965) * 0.98 for the 
* " * 
value of Q- obtained from Eq. 49. This value of 0 n is used to 
Umax ^ 2max 
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check whether the distribution of V 9 along r 9 , obtained from Eq. 47 

with a chosen value of meets the continuity requirements. To 

this end the ratios V thus determined, are introduced into Eq. 52 

and integrated from R ou /R to R 0 /R . If the resulting quantity 

^2max < ^ i ^^ ers f rom Q^ max /Rg of Eq* 54, the process must be repeated 

with a new value of V * until agreement is reached. However, the 

m2 * 

final distribution of along r^ must also produce the same radial 

shift of the streamlines as was determined from the first approxima- 
tion, and which was used to establish the function f(r 9 ) of Eq . 46 
* z 
and V 9 by Eq. 47. The equation of continuity is again used to obtain 

r 2 as function of r^ from the new ~ distribution. If this re- 

lationship between r^ and r^ differs from the one obtained with the 
data of the first approximation, it must be used to establish a new 
function f(r 9 ) of Eq. 46 which has to be introduced into Eq. 47 to 

calculate V 0 . This process requires the same iterations for V ^ 
a 2 m2 

as discussed above. The iterative procedure must be repeated until the 
functions r ^ vs. r^ assumed for the determination of f(r 2 ) and that 
obtained from the resulting profile of V ^ are identical. 

7. RESULTS OF FLOW CALCULATIONS 

* * 

Table V lists the ratios V * and V „ of Table IV, and additional 

al a2 9 

values at other radius ratios, which were used to calculate the 
* * 

quantities and in accordance with Eqs. 49 and 52, by integra- 
tions with the trapezoidal rule. It can be noted that the initial 

* 

assumption * V £ * does not satisfy continuity since 

(0.9832)Q. whereas by Eq. 54 with k = 0.98 it should be (1.0204) 

^ lmax y d 

^lmax ° r e 1 ua ^ to 0*404947. 

The ratios Q]/Qi max > anc * ^2^2max ^ a ^-*- e v are 1° 

Fig. 9. For the condition Q 0 /Q 0 * Q 1 /Q 1 , the radius ratio r 0 

2 2 max 1 lmax z 

for a stream surface that has the radius ratio r^ ahead of the rotor 
can be read on the upper horizontal scale of the figure. The relation- 
ship r^ vs. r^ thus obtained is shown in Fig. 10, which also gives a 
curve Ar * (r^ - r^) as function of r ^ in a larger scale drawn through 
the data points. From this curve were obtained the derivatives 
3(Ar)/9r2 ^ graphical means. Then 
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a(r 2 + Ar) 



• 1 + 



( 55 ) 




3(Ar) 

3r„ 



The numerical values of these derivatives and of r^ for different 
ratios r 2 are listed in Table VI, together with the function f(r 2 ) 
of Eq. 46 for X - 0.36578. 

The function f(r 2 ) is plotted in Fig. 11. Interpolated values 
between the calculated points were used to integrate f(r 2 ) with Simp- 
son's rule, i.e. by placing a parabola through three neighboring 
points. The results of the integration in accordance with Eq. 47 are 

listed in Table VI and plotted in Fig. 11. Table VI also shows the 

* 

calculated values of V 9 In accordance with Eq, 47, for several assumed 

values of V ^ > and the corresponding magnitudes of of Eq. 52 

obtained by integrations with the trapezoidal rule. As pointed out 

* * 
earlier, the value of at * 1.25, which is denoted by Q 2 max> 

should equal 0.404947. For V 0 * 0.838 this condition is closely 

m2 * 

realized and the corresponding values of V . are taken to satisfy 

* aZ 

the continuity requirement. For V j * 0.838 the ratios Q2'^2max are 
shown in Table VI also. 

It must now be checked whether the relation between r^ and r ^ 
from the first approximation of Fig. 9, which was used to calculate 
the velocity distributions of Table VI, agrees with the function r^ vs. r^ 
that is obtained from the second approximation. Since the change of 

it it 

V^ along r^, hence of V ^ along r^, is the same for both approxima- 
tions, this check can be made by plotting Q2^2max as ^ unction r 2 ^ or 
both sets of data. Figure 12 shows a curve of Q2^2max vs# r 2 < * rawn 
through the values of Fig. 9 which are marked by crosses. The data 

points given by circles in Fig. 12 represent the values of Table VI 
* 

for V^ * 0.838. It can be noted that these points lie very closely 
on the curve of the first approximation, hence additional iterations 
are unnecessary, and the data inside the heavily framed portion of 
Table VI can be considered to represent the actual flow conditions 
after the rotor for the set of initial assumptions; that is, if curva- 
ture effects and entropy gradients are ignored. This final distribu- 
* 

tion of V 0 is also shown by the dash-dotted curve in Fig. 8. 
a 2 
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Table VII summarizes the flow properties ahead of and after the 
rotor in accordance with the data of Table VI. Since the velocity 
diagram of the blading is that of Fig. 3, there are, with Eq . 41 



and 



Thus 



and 





ujR 

m 




(56) 





ujR 

m 




wR 

m 




(57) 



tanO^ 




(58) 






W 



ul 



V * 
al 



(59) 



With Eq. 45 and from Fig. 3 



and 



W 



u2 



!u2 

u>R 

m 



“ R 2 - V u2 

a) R 
m 



r - V 
2 u2 



tan( 



u2 



a2 



W 



tang. 



u2 



Also 



a2 



wR 



m 



al 

cosa. 



W, 



3 

k)R 



m 



al 

cosB, 



(60) 



(61) 



(62) 



(63) 



(64) 
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(65) 



V. 



U)R 



m 



a2 



cosa. 



W. 



s 

caR 



m 



a2 



cosB, 



( 66 ) 



Table VII also lists the NASA diffusion factors of Eqs. 17 and 

18, where the subscripts o in (D n ) and (D_) indicate that these 

K o oO 

factors are for a solidity a * a * 1 at R/R * 1, and that the chord 

m m 

of the rotor and stator blades does not change in radial direction. 

In actuality the diffusion factors should be calculated for the 
conditions along the stream surfaces by taking account of their 
radial shift. Then, by Ref, 3, for the rotor 



D 



R 




R 1 W U1 - R 2 W u2 
0^(1^ + R 0 )W 



2 1 



(67) 



A similar expression is defined in Ref. 3 for stators. However, 
in the present case where the radial shift is slight, the relations 
of Eqs. 17 and 18 can be used with very small error. In fact, the 
differences are so insignificant that the blade profiles in stator 
and rotor will be determined for cylindrical stream surfaces, with 
the actual flow angles of Table VII, instead of along the actual 
stream surfaces. This procedure avoids the difficulties associated 
with the definition of the blade shapes for manufacturing purposes. 

If the profiles are considered to lie on arbitrary stream surfaces, 
the resulting blade shapes must later be intersected by cylinders, 
or by planes tangent to them, to be able to obtain profile coordi- 
nates that can be used and checked during manufacture with reasonable 
cost. The available calculating methods for the profile shapes and 
their angular orientation are not accurate enough to warrant the 
difficulties associated with the afore-mentioned procedure, particu- 
larly not because the theoretical, axisymme trie , stream surfaces are 
very nearly cylindrical in the present case. It must also be recog- 
nized that the actual stream surfaces in a blade row are not axisym- 
metric, hence particles that move along the two sides of a blade 
travel along different paths with different radii so that the flow 
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does not have the same two-dimensional character that exists in 
rectilinear cascade test rigs which were used to establish the profile 
design data. 

In addition to, or as a replacement of the diffusion factor, the 
so-called equivalent diffusion ratio D is frequently used to estab- 
lish the blade loading limits and the profile losses, especially if 
the flow incidence angles differ from those for minimum losses. By 
Ref. 4, for rotors at the design incidence angle, and for flows with 
varying through-flow components on cylindrical stream surfaces. 



EQR 



w, l 1 - 



12 + 



0.61 



2 

cos 6 



2AW n 

u 

1 v . . v - 
al «f a2 , 



( 68 ) 



For stators at the same conditions. 



EQS 



I* r 

v i l 



1.12 + 



0.61 2 
cos a 



2AV 



2 V V 
al + a2 



(69) 



Experience has shown that the quantities should not exceed values 
of about 1.7 since incipient flow separation is likely to occur for 
D between 2.0 and 2.2. The subscripts o in (D_ rt _) and (D_ rt _) 

of Table VII indicate that a„ * o. * 1 at the mean radius R a nd that 

R S m 

the blades have constant chord at all radii. 

8. BLADE PROFILES AND STAGGER ANGLES 

The thickness distribution of Fig. 13 will be used for the 
profiles of the stator and rotor blades. It corresponds to a C.4 base 
profile (See Ref. 5), except for the portions near the entrance which 
were thinned to arrange a smaller leading edge radius. The camber 
line of the profile is supposed to be a circular arc. Figure 14 shows 
the method that will be adopted to determine the profile coordinates 
for known camber angles a and maximum thickness ratios t/c. The 
orientation of the blade profiles in the cascade is fixed by the stagger 
angle y as indicated in Fig. 14. 
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The profile data are determined with the two-dimensional method 
of chapter VI of Ref. 6, without taking account of the corrections 
for three-dimensional effects that are given in chapter VII of Ref. 6. 
It is felt that these corrections are not sufficiently well supported 
by experiments. Moreover, since the blades in the test compressor 
can be set at different angles the twist of the profiles along the 
blade height is more important than the stagger angle per se. 

The camber angle <j> for a rotor profile is obtained from 



* * 



Ag - i +6 
o o 

1 - m + n 



(70) 



where 



Ag - e x - b 2 



i 

o 



( Vsh <Vt 



( VlO 



(71) 



6 

o 



<V S h <K S>t 



<VlO 



(72) 



(i )= zero camber incidence angle for t/c ■ 0.1 

(6 q )^q ■ zero camber deviation angle for t/c * 0.1 

The recommended value of the correction factors (K. ) , and (K.) , is 

1 sn o sn 

1.1 for C.4 profiles. The design incidence angle i for minimum loss is 

i = i + n<j> (73) 

o 

and the deviation angle 6 for these conditions is 

6*6 + m<i> (74) 

o T 

The quantities (i ) , (6 ) , m, and n are functions of the inlet 

0 10 0 10 

flow angle B. and the blade solidity a. The correction factor (K ) 

1 t 

and (K^) depend on the thickness ratio t/c only. The blade stagger 
angle y of Fig. 14 is equal to 



Y “ 





(75) 
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The stator blade profiles are obtained with Eqs. 70 to 7 A also if 
8-^ is replaced by o^, and Ag by Aa ■ 02 - a^. The blade stagger angle 
is then 

Y - a 2 - i - £ (76) 

For particular conditions the method of Ref. 6 gives excessively 
large negative incidence angles. Although negative incidence angles 
are beneficial for good stall margins, they may cause increased losses 
at the design point because of separations on the pressure side near 
the leading edge, especially if the latter has a small radius. If 
the above-mentioned situation occurs, reduced negative incidence angle 
i* will be assumed. Reference 6 shows that for these conditions the 
corresponding deviation angle 6 is obtained from 

6’ - 6 + (i» - i) (||) (77) 



where the function (d6/di) depends on solidity and inlet flow angle. 

For a stator blade profile, for instance, there is 

Aa®4+i~ 6 * ^ + i* - 6* 

where the first expression on the right-hand side holds for the design 
conditions obtained with the usual approach, whereas the second 
holds for an assumed incidence angle i* from which the necessary blade 
camber angle <J>* can be determined. With Eq. 77 

> > f > > / JX \ 

<(> = Aa - i + 6 - Aa-i + 6 +(i - i) \"^±) (78) 

* 

For this modified camber, the blade stagger angle y of a stator 

is y * a- - i‘ - (79) 

2 

The axial distances between the blade rows of the compressor 
are 2.875 inches. The axial clearance between the rows must be about 
0.5 in. to provide room for flow survey probes, so that the axial blade 
widths must be less than about 2. A inches. Further, since each rotor 
has 30 and each stator 32 blades, the solidites of the rotor and the 
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stator cascades must remain within certain limits. 

The diffusion factors D and the equivalent diffusion ratios D 

LQ 

of Table VII are measures for the aerodynamic loading of the blade 
profiles in a cascade. Excessive loading is associated with flow 
separations, increased losses, and reduced efficiencies. At low 
blade loadings the efficiency of the cascade is decreased also because 
of large boundary layer losses. Optimum efficiencies of a row of blades 
are obtained if at all radii the blade profiles have the highest possible 
loading and the same margin of safety against stalling. The choice of 
the latter depends on the operating characteristics of the system for 
which the compressor is used. 

Table VII shows that this condition is not realized in symmetri- 
cal bladings with blades of constant chord. The diffusion factor of 
the rotor is highest at the outer radius and very much below permissible 
limits at the hub radius. The opposite trend would give a better per- 
formance since it is known that for the same diffusion factor the losses 

are higher at the rotor tip than at the hub. The data for the stator 

in Table VII show that the profiles are lightly loaded at the outer 

radius and have increased loading at the hub radius. For stators with 
radial clearance gaps at the hub radius the loading should be reversed 
also to obtain optinum performance, since the stator losses near the hub 
radius are larger than those at other radii even if the diffusion factor 
were constant everywhere along the blade. 

Within the limitations imposed by the rotor and stator blade 

numbers and the necessity of maintaining sufficiently large axial blade 

clearances for the arrangement of flow survey probes, the blade chords 

of rotor and stator will be varied along the radius to slightly improve 

the distributions of the diffusion factors of Table VII. The chosen 

chords for the rotor blades are given in Table VIII, those for the 

stator blades in Table IX. From Table VIII it can be seen that, whereas 

for constant chord the factor D for the rotor changed from 0.2062 to 

K 

0.4376 from hub to tip, it now has values of 0.2498 and 0.4274 at these 
respective stations. The stator solidities obtained with the chosen 
chords are generally lower than those for constant chord in Table VII, 
which gave variations of from 0.4107 to 0.3326 from hub to tip. By Table 
IX these respective values are 0.4240 and 0.3762 for the stator blades 
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with varying chord. These data show that with the symmetrical blading 
of Fig, 3 the earlier-mentioned optimum blade loadings cannot be achieved 
by varying the blade chord. Symmetrical bladings produce favorable 
conditions as far as the uniformity of the inlet flow velocities along 
the radius is concerned. Table VII shows for instance that the relative 
rotor inlet velocity decreases from 100 percent at the hub to 95.5 
percent at the tip. The stator inlet velocity decrease in radial 
direction also, namely, from 100 percent at the hub to 81 percent at 
the outer radius. The actual magnitudes of and are nearly equal. 
Hence, in contrast to other blading types, symmetrical bladings are 
not Mach number limited by the conditions at the outer radius. For a 
given maximum Mach number, they can operate at higher peripheral speeds 
than other bladings and are therefore capable of producing a higher 
pressure rise per stage. 

The best design of a blading would however be one that has uniform 
and equal velocities and along the radius, which are small with 
respect to the peripheral speed, and where the diffusion factors for 
all profiles are as high as the necessary stall margins permit. 
Investigations whether such bladings are possible will be carried out 
with a new computer program that is based on the calculating approach of 
this report, but with the additional refinement that the effect of the 
different losses along the radius will be taken into account. 

The profile thicknesses of Table VIII and IX were chosen for practi- 
cal and aerodynamic reasons. For the rotor the blade thickness increases 
about parabolically from tip to hub to obtain reduced centrifugal 
stresses, and the thickness of the stator blades varies about linearly 
to obtain small ratios t/c at the hub radius, and sufficient blade 
thickness at the outer radius for the attachment of the blade to its 
root section with the circular platform that is shown in Fig. 9 of Ref.l. 

Whereas the incidence angles i of the rotor blade profiles of Table 

VIII seem reasonable, excessively large negative incidence angle are 

obtained with the design method of Ref. 6 for the stator profiles of 

Table IX at large radii. The data inside the heavily framed portion of 

Table IX are therefore replaced by those for the assumed incidence 
> 

angles i < The last lines in Tables VIII and IX give the axial widths 
of the rotor and stator blade profiles. For an axial spacing of 2.875 in. 
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of the stator and rotor rows, the average axial clearance between the 

rows is 2; 875- (1 +1 )/2. Hence this axial clearance is about 

axr axs 

0,47 in. at the hub radius and 0.76 in. at the outer radius. These 
clearances are large enough for blades with the usual radial stacking 
line through the centers of gravity of the different profiles. It must 
however be examined whether for blades with swept-back leading edges, 
which are to be tested during the program, the blade chords of Tables 
VIII and IX permit to arrange inter-stage instrumentation. 

The data of Tables VIII and IX will be used to establish the profile 
coordinates with the relations of Fig. 14. For each profile there must 
be determined its cross-sectional area, the coordinates of the center of 
gravity and the angles of the principal axes with respect to the refer- 
ence line of Fig. 14, and the minimum and maximum moments of inertia. 
These quantities are necessary to establish the stresses and the vibra- 
tory characteristics for particular three-dimensional arrangements of 
the profiles in the build-up of the blades. 

9 . COMPRESSOR PERFORMANCE . 



If the losses in the inlet guide vanes of the compressor are 
ignored, the total pressure P^ ahead of the first rotor at station (1) 

is everywhere constant. However the static pressure p^ at station (1) 
changes along the radius, since for incompressible flows 



tl 



- SL 
2 



(80) 



At station (2) after the rotor the total and static pressures are 

denoted by P^ and p^, 

where 



P - v 1 
t2 2 V 2 



(81) 



In contrast to prior usage in this report the region after the stator 
is called station (3). At this station exist the same velocities 
and flow angles as at station (1) but the pressures differ from those 
at (1) , and the static and total pressures p^ and P^ are related by 



> - £- v 2 
t 3 2 V 1 



( 82 ) 
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From Eq. 7(44) of Ref, 2 the conditions along a streamline in the 
rotor between stations (1) and (2) can be expressed by 



P 2 + 2 W 2 



P 2 



K - t> - , o tT ^ m 2 

‘ 2 “ R 2 * P 1 + 2 W 1 " 2 “ R 1 ‘ AP R 



(83) 



where AP^ represents the losses in the rotor. In Eq. 83 the differ- 
ences in the geopotential energies of the fluis are ignored. It is 
customary (See Ref. 6) to express the rotor losses by 

4P r ■ y rI w i 2 < S4 > 



In Ref. 3 the pressure loss coefficient Y is related to the diffusion 

K 

factor D and the rotor solidity a D by 
R R 



2g R 

cos6^ 



0.004 + 



0.0639(D + 

K 



2.91 

0 . 1 ) + 



2.02 



3.77 



0 . 228D t 



(1 -X r ) 



where 







( 86 ) 



The subscripts t and h refer to the outer and hub radius, respectively. 

Equation 85 is an approximate mathematical formulation of the empirical 

loss correlations of Ref. 6 which attempt to account for tip clearance 

effects with the factor X , which increases the losses toward the rotor 

R 

blade tip. 

In Ref. 6 the loss in total pressure in the stator blade row is defined 
by 



AP, 



P - P 
t2 t3 



Y — V 2 
*S 2 V 2 



(87) 



where by Ref. 3 the stator loss coefficient Y<, can be approximated by 



(85) 
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2 o. 



cosa 



1 



2.91 2.02 3.77 

0.004 + 0.0639 (D g + 0.1) + 0.057D g (1 - X g ) 



( 88 ) 



The quantities and are the stator solidity and the stator dif- 
fusion factor, respectively , and 



X 



s 




(89) 



Due to the ratio X^ the loss coefficients of Eq. 88 increase 
toward the hub radius where the stator blades have a radial clearance 
gap. Evidently Eqs. 84 and 88 do not take proper account of the tip 
clearance and end losses. A discussion in Section 13.7 of Ref. 2 
shows that these losses must depend on the blade loading, the ratio of 
tip clearance and blade height, and the ratio of blade spacing and blade 
height, in addition to solidity. The present research program endeavors 
to establish a better and more physically correct evaluation of the 
end losses and tries to investigate whether peculiar blade shapes can 
be used to minimise their detrimental effects on stage performance. 

From Fig. 3 
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83, and with Eq. 84, 
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With Eqs. 39, 40, 80, and 81 



P t2 - P tl * p “ 2 R « 2 X - Y R ! W l 2 



(90) 
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where X is a constant, equal to 0.36578, for the bladings under considera- 
tion. With Eq. 64 the above relations, rewritten, gives 
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From Eq. 87, with Eq. 65 
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and 



P p 

*t3 tl 

[ t3 2 _ 2 

pu R 

m 



Y 2 Y 2 

x - / (W x *) - ^ < V 2 * ) 



(93) 



For frictionless conditions the coefficients Y and Y_ would be zero 

R b 2 2 

and the theoretical total pressure rise in the stage would be Xpu Rm 
by Eq. 93. Hence the total-to-total stage efficiency becomes 



1 - 



V<>' 



+ vv 



2X 



(94) 



Equation 94 should actually be evaluated along particular stream 

* * 

surfaces; that is, and should be taken at the corresponding 

radii R^ and R 2 of a stream surface. However, because of the small 
radial shift of the flow paths in the stage, and for the other reasons 
explained in paragraph 8, Eq. 94 will be determined for cylindrical 
stream surfaces and constant radii R = R^ * R 2 . A volume flow averaged 
stage efficiency can be evaluated from the different efficiencies 
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n t along the radius by 
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( 95 ) 



With Eqs. 49 and 50, with r 1u ■ R 1U /R , r, . * R, /R , 

In J.n m it it m 
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* 

where from Table V, = 0.397007. The average total pressure 

at station 3 is then obtained from 
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In addition to the changes of the total pressures and * n 

radial direction, the distributions of the static pressures p^, p^> 
and p^ are of interest also for comparisons with the measured 
pressures and for the determination of the aerodynamic blade forces • 
From Eqs. 80 and 63 
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From Eqs. 81, 91, and 65 
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From Eqs. 82, 93, and 63 
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Table X shows the loss coefficients in rotor and stator for the 
solidities of Tables VIII and IX, and gives the dimensionless pressure 
coefficients. It can be noted that highest efficiency * 0,9422 
occurs at the mean radius, whereas the efficiencies at the tip and the 
hub are 0.8238 and 0.8857, respectively. The last line of Table X 
lists the integrands of Eq. 96. By trapezoidal integration the inte- 
gral of Eq . 96 is 0.363035, hence the average stage efficiency is 



» 0* 363035 
n t " 0.397007 



0.91AA 



Then, by Eq. 97, with X * 0.36578 



— t3 - tl 

n t 3 " m °- 33448 

pw R 

m 



Figure 15 is a plot of the pressure coefficients of Table X to 

show the pressure changes along the radius of the stage. It can be 

noted that equal increases of the static pressures in rotor and stator 

occur at a radius ratio r «* R/R of about 1.1. At smaller radii most 

m 

of the static pressure rise occurs in the stator, whereas the opposite 
occurs between r ■ 1.1 and the tip. Hence because of the resulting 
radial velocity distributions a symmetric blading does not have a 
constant degree of reaction of 50 percent. 
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With the established data the operating conditions of the compressor 
can be determined so that 135 HP are absorbed by one, two, or three 
stages. By Eqs. 30, 31, AO, and 50, the driving moment M for one stage 
is 



2 2 

M - Q. p a>R X - 2ir w R 
lmax o m m 



p Q. X 
o lmax 



where Q * - 0.397007 from Table V. Then, by Eq. 36, for 135 HP, 

lmax _ , 

R - 1.2 ft, and p * 2.371(10 ) slug/ft , 

m o 



u> 



(135) (550) 



13.793 (10°) 



2ir 



p R 
o m 



*1 max 



or 



(u ■ 239.82 radians/sec 

and N - u » 2290 rpm 

Then, by Eq. 50, 

Q, * 103A ft^/sec 
X lmax 

From Eq. 97 with X ■ 0.36578 and ■ 0.911A 

P . - P - p o» 2 R 2 X ~ - 62.4 lb/ft 2 - 12.6 in. H.O 
1 3 tl o m t l 

The following table gives these data if one, two, or three stages 
absorb a driving power of 135 HP. 
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One Stage Two Stages Three Stages 







Z - 1 


Z = 2 


Z - 3 


Speed N 


(rpm) 


2290 


1818 


1588 


Volume Flow Rate Q, 

lmax 


(f t^/sec) 


1034 


820 


718 


Pressure Rise ZCP^-P^) 


(in. H 2 0) 


12.6 


15.9 


18.1 


Rotor Tip Speed ojR^ 


(ft/sec) 


359.7 


285.5 


249.4 



Before deciding on a drive for all three speeds, a study will be 
undertaken to determine the blade and rotor stresses, and to investi- 
gate whether critical shaft speeds and bearing capacities permit 
trouble-free operation for all cases. 
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FIG. 2 VELOCITY DIAGRAM OF SYMMETRICAL 
STAGE FOR FLOW ON CYLINDRICAL 
STREAM SURFACE WITH R - R 
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FIG. 



VELOCITY DIAGRAM OF STAGE WITH 
FLOW ON ARBITRARY AX I SYMMETRIC 
STREAM SURFACE 
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FIG. 5 VOLUME FLOW RATE Q OF COMPRESSOR FOR POWER OF 135 HP 
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FIG. ft SPEED OF ROTATION N OF COMPRESSOR FOR POWER OF 135 HP 



SYMMETRICAL ft LAD IMG WITH • - I 





FOR ONE STAGE 



FOR TWO STAGES 



FOR THREE STAGES 



4 




FIG. 7 FLOW DEFLECTIONS IN ROTOR CASCADE AT DIFFERENT VALUES 

OF FLOW INLET ANGLE a AND NASA DIFFUSION FACTOR D 

1 rn m 

AT BLADE SOLIDITY a - 1 

tn 







film 

37 






M.W.V. 



i 




lucivi vgurtt fLow 
! FCRST1 AmcetlHMIi 



AHEAD OP 



. ittfa! Rdi lap 



-Fttt' 



FIG 



head 



e a 



FerC4F.C«g« of Tokel Vokumfc 



^ 2 m»x 



< 3 . 



r -2 



b«t**etL Ladiua E 



and 



it 



man 



* 06 ! Mti 0 .^ 12 , 



1,4 



t.aw 



1 M 



± 5 , 



0 



‘ma 






*~r *~r 



*-H- 



-►«-r 









5 



39 



M. MV. SfzofTo 





tfiKii 



AO 








M.H V. 8/2i/lo 



FIG. 14 DETERMINATION OF PROFILE COORDINATES 




GIVEM: CHORD Cf MAX. THicKMESS k; CAMBER ANGLE <f; STAGGER AM6LE ^ 

c 



R - — 

c_ 2 sinCT/O 1 



Sin © 



~*r~ •('-¥■) B) 



zt 



* = _L 

C IOO 



. y c 7 t 

* c IOO c 



xl* X + V sin 0 5 c X ?_Jf 1 si n C ^A ) 

J \oo ioo c ' loo / 

u. v 




T X 



l }■ ifc •-"*«>] -rtscw H ' c “ 



y</ c ; 

y./c J 


|. cos ©| 


[ 1 7 7 il 

l 2 sin (•%.) + 100 c J 


- icotfy*; 



M M.V. 8/29/70 



44 



FIG. 15 CALCULATED PRESSURE DISTRIBUTIONS IN STAGE (SEE TABLE X) 

P t - Total Pressure , p - Static Pressure 
R - Radius, t» - Angular Velocity, p . Mass Density 
Subscripts: (l) Ahead of Rotor 

(2) After Rotor, Ahead of Stator 

(3) After Stator 

(m) Arithmetic Meah Radius of Blading 
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TABLE I 



THREE-DIMENSIONAL FLOW CONDITIONS IN COMPRESSOR STAGE WITH VELOCITY DIAGRAM 
OF Fib. 2 d)N CYLINDRICAL STREAM SURFACES 



DATA AT MEAN RADIUS R 



m 



Am = 
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m 


m 
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UJ 
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r^ 
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TABLE I 

THREE-DIMENSIONAL FLOW CONDITIONS IN COMPRESSOR STAGE WITH VELOCITY DIAGRAM 
OF FIG. 2 ON CYLINDRICAL STREAM SURFACES 



DATA AT MEAN RADIUS R 



m 



0 



im 



4-S* 



O. 35 
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CT l o 



Y ss . 4-1 4-fo2.(e<t. i») 
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IX 
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TABLE ic 

THREE-DIMENSIONAL FLOW CONDITIONS IN COMPRESSOR STAGE WITH VELOCITY DIAGRAM 
OF FIG. 2 ON CYLINDRICAL STREAM SURFACES 



DATA AT MEAN RADIUS R 

IT) 




40* 



45RII 



D = o -4o 







- I ° 

. 755^4 



ta r 




0.75 


O. 075 


1 . O 


f . 115 


1 . 25 


\ o 




.41743 


.4© |ol 


,84944 


.81191 


•7^464 


•X 


Vaa /? w ) 


1-05 2o3 




.*6444 


.7S4>6 


.6 loll 


13 


trfn y 3 , 


.14146 


.77 674 


.637 1 » 


93150 


1.06444 


14 


t-r» 


ot 552 




.31 106 


•47577 


.723 37 


16 


k 43 n ^ 2 , 


.6*737 


.71 » 99 




1 . 01 6 1 6 


1-325 36 


1 5 


too oC) 




.154 39 




.435 SB 


58 097 




fit 


36*35' 


37*5o’ 


40* 


42* 58 ' 


46° 47 ' 




/** 


3 c 45* 


10 * 16 * 


17 * 16 ' 


25*25’ 


35* 53’ 




cL x 


32*55' 


35*54’ 




45*28' 


52*58 ’ 




4?C| 


4°lft' 


11 * ©' 




23* 32' 


3o*o3 * 




CO 6 /3, 


.So 3 oO v 


*78 975 


.76604 


.73172 


.48476 






.44785\ 


.4B398 


. 454-57 


.40322 


.6)023 




CoS <*1 


.63445 . 


.©' ©#o 




.70142 


,60230 




CoS 


.44 719 


•38165 




.41680 


.86467 


Zl 




.2704-5 


.340*1 


.4 


.46148 


. 5282* 


22 




. 44 066 


.42 325 


. 4 


. 37 85 4 


.35728 




IV, /(to R„) 


l.l47>o| 


f.l4o66 


U349B 


1. 12464 


J.10452 




W a /w, 


.41236 


.86 109 


. 80224 


. 7426 3 


.67864 




V*/(to R„) 


1.25324 


103-541 


1.13*70 


l.o7 547 


1.0) 24 8 




v,/v. 


.73 45 1 


.76910 


.©0224 


.8 346 1 


.86733 






32 0 So 1 


27*34’ 


2 2 *44' 


1 7 *lV 


10 * S4 1 




| 

Act - _ <A ( 


26 °37' 


24*54' 


22* 44' 


2 1 °S6' 


12 *49' 



48 



M.M.V. fi/23/lo 



table 12 



THREE- DIMENSIONAL FLOW CONDITIONS IN COMPRESSOR STAGE WITH VELOCITY DIAGRAM 
OF FIG. 2 ON CYLINDRICAL STREAM SURFACES 
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TABLE 1 FLOW RATE CALCULATIONS FROM FIRST APPROXIMATION 
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TABLE YDI RCdOR BLADE PROFILES FOR ASSUMED CHORDS AND THICKNESSES 
NUMBER OF ROI\)R BLADES - 30 
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TABLE X ROTOR AND STATON LOSSES , PRESSURES AND STAGE EFFICIENCIES 
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with symmetrical bladings for a research program to investigate tip clearance 
effects in the three- stage compressor of the Turbo- Propulsion Laboratory, NPS. 

It establishes the blading data and the stage performance with an iterative 
three-dimensional approach, and gives design criteria for the drive and the flow 
measuring device of the test unit. The calculated distributions of the flow 
properties in the stage will be used for future comparisons with test data. 
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